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We study tunneling in semiconductor heterostructures where the constituent 
materials can have a direct or indirect bandgap. In order to have a good de- 
scription of the lowest conduction band, we have used the nearest-neighbour 
sp 3 s* tight -binding model put forward by P. Vogl et al.. A recursive Green- 
function method yields transmission coefficients from which an expression 
for the current density may be written down. The method is applied to 
GaAs/AlAs heterostructures. Electrons may traverse the AlAs barriers via 
different tunneling states ipr an d ipx (TX mixing). With an applied bias 
V^0.5 V electrons may enter the GaAs collector contact in both the T and 
the X valley (TX transfer). We have studied a number of GaAs/ AlAs struc- 
tures. For very narrow barriers there is little TX transfer, but AlAs barriers 
wider than about 25 A act as 'TX filters", i.e., most transmitted electrons 
have been transfered to the X valley. 
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I INTRODUCTION 



Tunneling in semiconductor heterostructures has attracted considerable 
interest over the last decade. A very important motivation factor has been 
the progress in advanced crystal-growth techniques like molecular-beam 
epitaxy (MBE). The ability to grow nearly perfect layered structures has 
enabled experimental verification of predictions based on relatively simple 
theoretical models. There is also a great interest in making electronic devices 
based on such structures. 

Most treatments of transport through heterostructures have been based 
on effective-mass theory. This is a good approximation when the different 
materials which form the structure all belong to the same category with re- 
spect to type of energy bandgap, direct or indirect. An example of a system 
with only direct-bandgap constituents is GaAs/Al x Gai„ x As, provided the 
Al concentration is sufficiently low, x^SOA. In a typical experiment, doping 
in the GaAs contacts yields a Fermi level Ep of the order of 10-50 meV 
above the conduction-band minimum E^ at the V point of the Brillouin 
zone (k = 0). Since the conduction band is nearly parabolic in the range 
from to Ep, all the incoming electrons are well described by a single 
effective mass mp(GaAs). Furthermore the lowest tunneling barrier is de- 
termined by the conduction-band minimum in Al^Gai-^As, which is also 
at the r point. Thus the tunneling states, through which electrons with 
energy E can traverse the barrier region, are characterized by an imaginary 
wave vector k = in, where k is determined by the tunneling-barrier height 
E^ (ALjGai^As)— E and the effective mass mp(Al a .Gai_a;As). 

For Al concentrations x > 0.4, Al^Gai-^As becomes an indirect-bandgap 
material. The valence-band maximum is still at the T point, but the 
conduction-band minimum is now E* , close to the X point, k = ^(100), 
at the edge of the Brillouin zone. Here ai is the lattice constant of the 
zincblende material. When indirect-bandgap Al^Gai-^As is used as the 
barrier material, it is no longer sufficient to take into account only the tun- 
neling states which correspond to the conduction-band minimum at V. Im- 
portant are also the states that correspond to the analytical continuation 
into the energy gap of le&l-k states at the conduction-band minimum near 
the X point. 'T states" and ll X states" may have comparable decay lengths 
since the lower barrier height of the X states, E^ — E, is compensated 
for by a higher effective mass. The importance of X states in tunneling 
through indirect-bandgap barriers has already been appreciated in several 
experiments and it has also been studied theoretically by a number 
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of groups 

In the present paper we want to focus on the X states in both the barrier 
and the contact material. It is well known that electrons in bulk GaAs 
may be transferee! from the T to the X valley if they are accelerated in 
a sufficiently strong electric field. This is the Gunn effect which may be 
accompanied by a region of negative differential conductivity in the current- 
voltage curve due to the low mobility of electrons in the X valley. The U TX 
transfer" is usually stimulated by some kind of scattering mechanism, e.g., 
elastic-impurity scattering. Our aim is to show how a heterostructure with 
one or more indirect-bandgap barriers can be used to control the TX transfer 
in a material like GaAs. Let us assume that the transport takes place in 
the (100) direction. Then, with a bias V > [£^"(GaAs) - E F ]/e applied 
across the heterostructure, incoming electrons in the left contact can end up 
in two final states in the right contact, in the T valley or in the X valley. 
The two outcomes are characterized by transmission coefficients Tp and Tx , 
respectively, which depend on the energy of incoming electrons, applied bias, 
and barrier parameters. In Sec. V below we shall explore these dependencies 
for GaAs/AlAs heterostructures in order to find effective "TX filters", i.e., 
conditions under which Tx 3> Tp. 

In transport experiments the measured quantity is usually the electric 
current vs applied voltage. Based on simple arguments one can write down 
an expression for the current density in terms of transmission coefficients 
and a factor ensuring that initial states are occupied and final states empty. 
In the present case a sum over different final states at equal energy au- 
tomatically provides a decomposition of the current density in a "direct" 
component Jr and a "transfered" one Jx- This simple procedure is indeed 
applicable in the present study of perfectly layered heterostructures where 
we ignore effects of disorder and inelastic scattering. If such effects are to be 
included, one must resort to a more general approach. In an Appendix we 
derive an expression for the current density which provides the starting point 
for extensions to more realistic calculations [19|. The general expression in- 
volves nonequilibrium Green functions and gives the total current density J, 
not the decomposition into Jr and Jx- However, in the case where disorder 
and phonons can be ignored, we will show how the general expression for 
J{V) can be decomposed and proven to be identical to the result which was 
written down directly in terms of transmission probabilities. 

The transmitted electrons will be subject to elastic and inelastic scat- 
tering and eventually come to thermal equilibrium somewhere in the right 
contact. This process may be described by scattering rates ^rr> r xx and 
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T xri w here the former two account for intravalley scattering and the latter 
describes relaxation of X electrons to the T valley. Clearly, txv sets the 
time scale for observing or taking direct advantage of having electrons in 
the X valley. 

The rates of elastic and inelastic scattering depend mainly on impurity 
concentration and temperature, respectively. At low temperatures the dom- 
inating inelastic process in GaAs is spontaneous emission of LO phonons. 
A rough estimate, based on Refs. [20, |2^], yields an inelastic relaxation 
rate on the order of 10 13 s" 1 in both the T and the X valley. A similar 
estimate yields an intervalley scattering rate t^y ~ 10 12 s -1 , taking into ac- 
count that electrons in the X valley may relax to the T valley via emission 
of LO phonons or via scattering off charged impurities (assuming impurity 
concentrations of about 10 17 cm -3 ). This implies that, with an average drift 
velocity of about 10 5 m/s, a transfered electron travels typically a distance 
of 1000 A in the X valley. 

One way of detecting X electrons, then, could be by means of a magnetic 
field B applied parallel to the crystal-growth direction, and taking advantage 
of the difference in Landau-level splitting (AEll = = TieB /m*) in the 
r and the X valley due to the difference in effective mass (m x 3> mf). 
However, in order to test predictions for the TX transfer, it may even, for 
carefully designed heterostructures, be sufficient to measure the total current 
vs applied voltage. If the nonlinear structure in the measured total current 
is well described by the theory, it seems reasonable to assume that the 
calculated decomposition into Jr and Jx also agrees with reality. Finally, 
the TX transfer may be checked by investigating the noise spectrum of the 
current. Quite often the noise spectrum of a given physical quantity may 
reveal sharper features and more information than a measurement of the 
physical quantity itself [23|. 

The effects we want to study involve states far from the local minima 
of the conduction band and also states in more than one valley. In order 
to describe the band structure correctly, at least qualitatively, we apply 
the semiempirical sp 3 s* tight-binding (TB) Hamiltonian put forward by 
Vogl et al. f24j. A feature of this model, of particular interest here, is 
the correct description of the transition in Al x Gai_ x As from a direct to an 
indirect bandgap. In addition it is a three-dimensional model which makes 
it suitable for evaluation of current-voltage characteristics. The sp 3 s* model 
was also used by Cade et al. 0j and by Yamaguchi [ 25 1 . 

We have organized the paper as follows. In Sec. II we describe briefly 
the sp 3 s* TB Hamiltonian. In Sec. Ill we derive expressions for the scatter- 



4 



ing coefficients. An expression for the current density is presented in Sec. 
IV, and in Sec. V we perform a systematic study of the TX transfer in 
GaAs/AlAs heterostructures. We conclude in Sec. VI and discuss briefly 
possible extensions and other applications of the present model. 



II. THE sp 3 s* HAMILTONIAN 

The semiconductors that we want to describe have the crystal structure 
of zincblende. The valence bands of these materials are usually well described 
by an eight-band sp 3 nearest-neighbour semi-empirical TB Hamiltonian 
p6| , |24| , the eight bands arising from the inclusion of four atomic orbitals 
on each of the two types of atoms, cation and anion. However, for electronic 
transport and tunneling we need a model which describes the lowest con- 
duction bands adequately, and for that purpose the nearest-neighbour sp 3 



model generally fails. In particular, as shown by Chadi and Cohen [27], it 
cannot produce an indirect bandgap in materials like AlAs, and this is the 
single most important feature required for our purposes. 

Vogl et al. [24] have overcome this deficiency by the ad hoc inclusion 
of an excited s state on each atom. The main effect of coupling these s* 
states to nearest-neighbour p states is to repel the p-like conduction-band 
levels near the X and L points to lower energies, thus producing the desired 
indirect bandgap. The resulting ten-band sp 3 s* Hamiltonian has thirteen 
independent TB matrix elements which are determined by fitting band- 
structure data. 

The independent TB matrix elements are given in a basis of symmetri- 
cally orthogonalized atomic orbitals \nbR), also called Lowdin orbitals pqj . 
Here R is the position of the atom, n denotes the type of orbital (s, p or s*) 
and b the type of atom (c for cation and a for anion). In Ref. p|| the Hamil- 
tonian for a bulk material is expressed as a 10 x 10 matrix in a basis \nbk) 
which is obtained by a discrete Fourier transform of the localized orbitals 
\nbR). The systems that we want to study are translationally invariant in 
the "parallel" plane (y,z). However, the heterostructure breaks the invari- 
ance in the crystal-growth direction which is usually chosen to be along the 
(100) axis. Thus it is convenient to represent the Hamiltonian in a basis 
\nbjk\\) with parallel momentum ku = (k y ,k z ) and "perpendicular" position 
Xj = jai/4: as parameters. The distance between nearest-neighbour planes 
is one forth the lattice constant ai ■ A simple inverse Fourier transform (with 
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the same normalization convention as in Ref. [p4|]), 



\nbjk\\) = L 



-1/2 
BZ 



dkx g 



-ik x ja L /4 



nbk), 



(2.1) 



yields the desired basis. Here Lbz = 8n/aL is the length of the one- 
dimensional (ID) Brillouin zone over which the k x integral is taken. 
In the new basis the Hamiltonian takes the form 



H 



Vca E c 

u a , 



U ca 

E a 
Vca 



V ac 



U nr E„ V n , 



(2.2) 



where each element is a 5 x 5 matrix. The matrices E c and E a are di- 
agonal and represent the orbital energies on cations and anions, respec- 
tively, whereas the "hopping" terms U ca , U ac , V ca and V ac involve the various 
transfer-matrix elements between orbitals on neighbouring sites. Obviously 
U ac = U\ a and V ac = V^ a so that the Hamiltonian is hermitian. In Ap- 
pendix A the elements of H are given in detail, and we have also included 
the numerical values of the matrix elements for GaAs and AlAs, taken from 
Ref. W4 



The Hamiltonian in Eq.(2.2) is formally identical to that of a two- 
atomic ID chain with interatomic separation a^/A and periodicity a^/2, see 
Fig. 1. 

Alloys will not be treated explicitly in the present paper. We only men- 
tion that the simplest approach would be the virtual-crystal approximation 
(VCA). This means that the Hamiltonian of an alloy A^Bi-^C is approx- 
imated by the weighted average of the Hamiltonians of AC and BC, e.g., 
U ca [A x Bi_ x C]= xU ca [AC] + (l — x)U ca [BC]. Scattering due to disorder in 
the cation planes is neglected in the VCA. That effect can be included by 
replacing VCA with the socalled coherent -potential approximation (CPA). 

Another effect which is ignored in this work, is elastic scattering due 
to interface roughness. Like alloy disorder, interface roughness breaks the 
translational invariance in the parallel plane. Hence the parallel momentum 
is no longer conserved, and the simple treatment of the next section must be 
modified. Interface-roughness scattering can also be accounted for within 
the coherent-potential approximation. 
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Ill SCATTERING COEFFICIENTS 



When the physical system is described with an effectively ID nearest- 
neighbour TB Hamiltonian, the scattering states can be evaluated efficiently 
with a recursive Green-function technique. This method has been applied 
to various problems |2^]. The present treatment will essentially be a gener- 
alization of the simple one-band ID chain described in detail in Ref. |3(| . 

Since we will discuss only processes where the parallel momentum ku 
is conserved, we will suppress fen in most of the notation of this section. 
It is furthermore convenient to label a cation layer and its neighbouring 
anion layer to the right with a common "site" index j. Assume that the 
heterostructure and the undoped spacer layers, if any, are located on the 
sites from j = 1 to j = N. Then the homogeneous "leads" extend from 
j = — oo to j = and from j = N + 1 to j = oo. Write the total wave 
function |^) as a sum of two pieces, 

\V) = \S a ) + \(f>). (3.1) 

The first term is taken to be the incoming plane-wave part with wave number 
k a = (k a ,k^) and energy E, restricted to layers up to and including the 
cation part of site 0: 

\S a )= £ e lk ^ 2 \j,a). (3.2) 

j<0,c 

The local part \j,a) is a 10 x 1 vector with coefficients a n b (n = s,p x ,p y ,p z , s*; 
b = c,a). These coefficients are determined by the solution of the Schrodinger 
equation for the bulk material at energy E and wave vector k a . The problem 
is now reduced to finding the "remainder" \<p), and since {E — H)\^) = 0, 
we have 

\<t>) = G R {E)[-{E-H)\S a )]. (3.3) 

Here G R {E) = lim r ^ +(£' — H + ir])^ 1 is the retarded Green function 
which also depends on fey via the Hamiltonian. From the definition ( |3.2| ) 
of \S a ), and since H only couples nearest-neighbour atomic layers, the term 
in square brackets in ([O]) is clearly zero for sites j > 0. Furthermore the 
Schrodinger equation ensures that (E — H)\S a ) = for sites j < 0, and one 
is left with nonzero terms only on the cation and anion layer at site j = 0. 
We may finally write 

\4>) = —G R (E)v\0,a), (3.4) 
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where v resulting from ( |3.3j ) reads 



4h 



iU ca 
-iU ac 



(3.5) 



The operator v can be interpreted as a local velocity operator. This is 
consistent with the velocity operator obtained from the Heisenberg equation 
of motion, v = — ^[x,H], with the position operator 



E 



|jc)^-(j c | + |ja)(— + -j){3a\ 



(3.6) 



Here we have split the vector |j) into its cation and anion parts \j c ) and \j a ), 
respectively. 

It is now straightforward to determine the scattering amplitudes. The 
total wave function may be written as 



I*) 



E 

3<0 



Jk a ja L /2 



+ 



+ 



E E 

j>N+l (3 
E EXa«b> 

l<j<N K 



/3<Z 



-ik ja L /2 



P 



Va f Jkpja L /2 
V P 



(3.7) 



The incoming plane wave with wave number k a may be reflected or trans- 
mitted with wave number kg, with the restriction that the total energy and 
the parallel momentum are conserved. The velocity factors are included so 
that particle conservation yields a natural probabilistic interpretation of the 
transmission and reflection coefficients, 



E [iW*! 2 + l r «/3l 2 



1. 



(3.8) 



Information about the various tunneling states are contained in the coef- 
ficients |x^ K | 2 ) the j dependence of which yields the spatial-decay rate of 
the state with (in general) complex wave number k K . An expression for 
the transmission amplitude t a p is obtained by comparing the projection 
{N+l,p\V) ofEqs. (O) and " 



(3.9) 
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In a similar way the projection {0, f3\^) yields the reflection amplitude 



r a(3 




4ih 



{0,P\G R v\0,a) -5 aP {0 a ,a\0 a ,a) 



(3.10) 



where the second term comes from the fact that \S a ) does not include the 
anion layer of site j = 0. Eqs. ( |3.9[) and ( 3.1C| ) are readily seen to reduce 
to the expressions for t(k) and r(k) given in Eq. (2.9) of Ref. [ 30 1 for the 
simple case of a one-band ID model. 

Note that the "left" states {j, j3\ in (3.9) and ( 3.1 0| ) are not simply the 
hermitian conjugate (j, j3\ of \j,/3). The reason is that the local projec- 
tions | j, (3) on a site j do not correspond to an orthogonal basis. Hence, 
(jj 0) 3(3/3' j an d we must construct a new basis of left states such that 
{j>/3'\j,P) = Sgai. Let B denote the 10 x 10 matrix whose columns are the 
ten states |/3) = |j, /?), independent of j. Then the matrix inverse B^ 1 is 
precisely the orthogonal basis that we need, with rows {/3'\. In addition to 
orthogonality one has, by construction, the closure relations 



& (3 



(3.11) 



We can bring Eqs. Q and (|3~10| ) into more symmetric forms by using the 
identity 

v\a) = v a \a}. (3-12) 

Here v a is the velocity [in the (100) direction] associated with the state |a). 
Eq. (fnp is a result of (pl^ ) in Appendix B. Using ( |37T^ ) in (|^) and 
(|3l0|) yields 



J-a.f3 
R a f3 



\ta(3\ 



^^^{iV + l,/?|G R |0,a}p 



\r a p\ 2 = \&y/\vp\v a {0,P\G R \0,a} - 5 a p{Q a , a\0 a , a}\ 2 . 



(3.13) 



IV. ELECTRIC CURRENT 

Transport in solids can be treated on various levels, depending on the 
approximations that are assumed at the outset. Examples are the semiclas- 
sical Boltzmann equation, the Kubo formula for linear response, and the 
Landauer formula which expresses the conductance of a system in terms of 
the single-particle transmission coefficients pi]]. 
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In the case of a perfectly layered heterostructure it is well known that an 
expression for the electric current can be written down directly in terms of 
the transmission coefficients. Assume that one can define chemical poten- 
tials [II and hr in the left and right leads, respectively. Under equilibrium 
conditions fj,i = fj,R, and when a bias V is applied between the two leads, one 
has /j,r = hl — eV. Strictly speaking, these potentials do not describe the 
situation close to the barrier structure, a region which is out of equilibrium 
because of tunneling electrons. Thus, [1l and /j>r correspond to asymptotic 
distributions in the leads. In the left lead, states with energy E are occupied 
with probability Jfd(E — fj,^) = [exp((E — h^/IzbT) + l] -1 . Here fpD is 
the Fermi-Dirac distribution function, ks is Boltzmann's constant, and T 
is the temperature. An electron in initial state \a) with velocity v a in the 
(100) direction has probability T a p of being transmitted into the right lead 
in a final state \/3), which has probability 1 — Jfd(E — hl + eV) of being 
empty. Analogous arguments hold for the states in the right lead. In order 
to find the total current flowing in the (100) direction, we must integrate 
over all possible states, i.e., over the Brillouin zone. In addition we have to 
sum over the different energy bands for each k. 

In the present case we have scattering processes that conserve parallel 
momentum and total energy, whereas the momentum hk in the (100) direc- 
tion may be altered. It is therefore convenient to replace the integral over k 
and the sum over energy bands by an integral over energy and a sum over 
momenta hk a : 

^ f dk ^ fdE fdEy 1 
Since {dE/dka)^ 1 = (7w a ) -1 , the total electric-current density is given by 
J = -{Y.j ^Yij ^T^E^ifME-^-fFDiE-VL + eV)}. 

a/3 ^ ' 

(4-2) 

Here e is the absolute value of the electron charge so that V > results in a 
net flow of electrons from left to right, i.e., a negative electric current. The 
sum over a and j3 is restricted to states with velocities v a ,vp > in the 
(100) direction, and we have taken advantage of time-reversal symmetry 
under which T(k a — ► kp) = T(—k/3 — ► —k a ). The integral over parallel 
momentum fey runs over the two-dimensional Brillouin zone which for the 
[100] plane in zincblende structures is a square with corners (0, ±27t/ol) and 
(±27r/a L ,0). 
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Note that the simple arguments above are no longer applicable if one 
wants to include effects of alloy scattering or electron-phonon interactions. 
In that case a full Green-function treatment is required, and in Appendix 
B such an approach is presented. For the ideal case, where (|4.2| ) is valid, 
we demonstrate explicitly the equivalence between the Green-function ap- 
proach and the transmission-probability approach. 



V. Results and discussion 

The model and method described above can be applied to study elec- 
tronic transport through arbitrary heterostructures with zincblende-type 
constituents. In this section we shall concentrate on the TX mixing and 
transfer, and results will be presented exclusively for GaAs/AlAs heterostruc- 
tures. We have used a simple, linear relation, 

AEc(x) = 0.9x (eV) , (5.1) 

for the r point conduction-band offset between Al x Gai_ x As and GaAs. This 
is a compromise between various empirical relations quoted in the literature 
| p2| . In combination with the TB parameters of Ref. [£4| (see Appendix A), 
Eq. (5T) yields a discontinuity of about 0.16 eV between £",5 ( AlAs) , i.e., 
the X point conduction-band minimum in AlAs, and ££(GaAs). This is m 
agreement with recent experiments 0—0. 

We believe it is instructive to base the following discussion on parts of the 
complex band structure of GaAs and AlAs. In Fig. 2 the solid lines denote 
the lowest conduction band in the (100) direction, i.e., for ku = 0. For AlAs 
we have also plotted the two branches of complex wave vector corresponding 
to the analytical continuation below the conduction-band minima at the T 
and X points of the Brillouin zone . These branches are found by solving 
the Schrodinger equation for fixed energy E. In the present model there are 
10 solutions k(E) = k^ + in that make up the full complex band structure. 
In Fig. 2 we have only plotted the imaginary parts, ky and kx, of the two 
branches that are relevant to tunneling in a GaAs/ AlAs heterostructure. 
For an incoming electron at energy E in the T valley of GaAs, nr(E) and 
kx(E) determine the exponential decay of the tunneling wave functions Vt 
and ipx through the AlAs barrier (s). Clearly, the relative importance to the 
tunneling process of t/t and ipx will depend on which has the slowest decay 
in the barrier region. For energies below, but close to the conduction-band 
minima we have k v {E) ~ [2mf (E£ - E)/h 2 ] 1/2 and K X (E) ~ [2m x (E* - 
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E)/^] 1 / 2 , where m T and m* x are the effective masses in the T and the 
X vailey, respectively. Thus, although the "A barrier" — E is much 
lower then the 'T barrier" — E, the large value of m* x yields kx > 
for energies up to about 20 meV above the conduction-band minimum in 
GaAs. Furthermore it is not only the exponential decay of the tunneling 
wave functions that determine their contribution to transmission through a 
barrier. At the material interfaces plane waves in GaAs must be matched 
to the wave functions ipr and ipx in the AlAs barrier. From symmetry 
arguments we expect, for an incoming V electron, a larger matrix element 
for the matching to ipr than to ipx- These effects are illustrated in Fig. 3 
where we have plotted the absolute value of the tunneling states as a function 
of position for a GaAs/AlAs potential step. The solid lines represent states 
at an energy 10 meV above E v c (GaAs). The stronger coupling to ipr is 
clearly observed by the fact that |^r| S> \ipx\ in the first layer of AlAs 
(i.e., for Nb = 1). Furthermore, since Kx~^r for this energy, ipx decays 
slightly faster than ipr, and tunneling will predominantly happen via ipr- 
At an energy of 100 meV above E^ (GaAs) (dashed lines) we see again the 
stronger coupling to ipr at the material interface. However, in this case kx is 
considerably smaller than Kr, and beyond the 10th AlAs monolayer we have 
\ipx\ > |Vt|- Thus, depending on the width of the AlAs barrier, tunneling 
at this energy may take place via ipr, ipx or both. 

In a typical experiment the Fermi level is about 10-50 meV above £f(GaAs). 
With the aid of Fig. 4 we will now sketch a qualitative picture of the tunnel- 
ing of an electron through a single AlAs barrier, where a bias V is applied 
such that transfer to the GaAs X valley is possible. Let us assume that 
the energy of the electron is far from any resonances in the system. At 
the left interface (x\) the matching of the incoming plane wave to tunnel- 
ing states ipr and ipx may be described with matrix elements Mrr(xi) and 
Mrx(x\), respectively. We saw in Fig. 3, in accordance with expectations, 
that Mrr S> Mrx- Since there is a bias across the barrier, the decay pa- 
rameters Kr and kx will now be functions of position within the barrier. 
Between x\ and %2 (typically a very short distance), we have kx > k>t, 
so up to this point tunneling via ipr dominates strongly. However, beyond 
X2 we have kx < Kr, and the amplitude for tunneling via ipx is "catching 
up" relative to that of ipr- Between X3 and X4 the energy of the tunneling 
electron lies above (AlAs). This means that kx = 0, and the amplitude 
of ipx is not further reduced whereas ipr continues to decay exponentially. 
Finally, at the right interface (x^) matching of the two tunneling states to 
transmitted states in the collector contact may be described by matrix el- 
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ements M aj g (a?4), with a,f3 = T or X. As above we expect the coupling 
between states of different symmetry (M^x and Mxr) to be much weaker 
than that between "similar" states (Mrr and Mxx)- 

Based on this simple picture we can make predictions concerning the 
two contributions, Jr and Jx, to the total current density in a GaAs/AlAs 
heterostructure. First, there is no FX transfer for biases V < [E* (G&As) — 
Ep]/e. However, the direct component Jr may still be strongly affected 
by the presence of the tunneling state ipx, hi particular due to resonant 
tunneling via confined X states in the barrier. We will briefly come back 
to this below. For biases V > [-^(GaAs) — Ep]/e, TX transfer is possible, 
and the total current density is given by the sum of Jr and Jx- Their 
relative contribution is determined by two competing factors. On one hand, 
Jr is favored by the stronger coupling to ■i/'r at the barrier interface. On the 
other hand, Jx is favored by the slower exponential decay of ipx through 
the barrier. For very thin barriers the interface coupling is the decisive 
factor, and we expect Jr > Jx- With increasing thickness of the barrier the 
difference in decay rate will eventually yield Jx > Jr, and the crossover to 
a a TX filter" takes place around a barrier width L-px that will be estimated 
below via numerical examples. 

The qualitative discussion above can only be expected to be valid in the 
off-resonant tunneling regime. The resonant features of tunneling through 
indirect-bandgap heterostructures have been discussed extensively in the 
literature |L7], 18], and they have been observed very clearly in experiments 
[^, [Tel l. In general there are two types of resonant states in a GaAs/AlAs 
structure. With two or more AlAs barriers there will be one or more GaAs 
wells in which one finds the "normal" quasibound states er confined by the 
r profile of the conduction band. On the other hand, the X profile of the 
conduction band reverses the role of barrier and well material, and as a result 
there will be resonant levels ex confined to the AlAs layers. These states are 
degenerate with the continuum of the GaAs contacts and hence similar to 
the socalled Fano resonances in atomic physics [j34f| . Predictions concerning 
resonance positions and linewidths are to some extent model dependent, as 
discussed in Ref. [18]. Here we shall not elaborate further on the resonant 
structure of the current, although tunneling via resonant levels is indeed 
taking place in the structures studied below (only ex in the single-barrier 
case, both er and ex in the double-barrier case). 

In order to determine the conditions for having a l TX filter", we have 
calculated current-voltage curves for a number of single- and double-barrier 
structures. Explicit results will be shown for single barriers only. We have 
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used a Fermi level of 10 meV above the GaAs conduction-band edge. Fur- 
thermore the potential profile is taken to be flat throughout the GaAs con- 
tacts, with a linear voltage drop across the barrier structure. In Fig. 5 we 
have plotted the two current-density components Jr (thick solid line) and 
Jx (thick dashed line) for a single AlAs barrier of width 40 A, for biases in 
the range where TX transfer is possible [35], and at zero temperature. We 
have also included the corresponding transmission coefficients Tr (thin solid 
line) and Tx (thin dashed line) at the Fermi level and with fen = 0. This 
shows that, although Jr and Jx are given by integrals over E and fen of 
Tr and Tx, respectively, their dependence on applied bias is to a very good 
approximation reflected in a single value of the transmission coefficient. Of 
course, the integrals over E and fey tend to give a smoother behaviour of 
J(V) than that of T(V). In addition, resonance positions and relative peak 
values may be slightly shifted in J(V) when compared with T(V). However, 
since we are not interested in details in the current -voltage curves, the fol- 
lowing discussion may be based on the behaviour of Tp and Tx at E = Ep 
and fen = 0. 

From Fig. 5 we see that a single AlAs barrier acts like a good TX filter 
for a width of 40 A. As a measure of the TX transfer we have calculated the 
ratio 

fdVT x (V;E = E F ;k ll =0) 

Rrx = =r , (5.2) 

JdVT r (V;E = E F ;k» =0) 



with limits of integration at 0.5 and 0.8 V ||36|| . The result is shown in Fig. 
6 where -Rrx is plotted versus the number of monolayers of AlAs in the 
barrier. As expected there is only negligible TX transfer for very narrow 
barriers, but already at 10 monolayers there is almost an order of magnitude 
more electrons being transmitted in the X valley than in the T valley. Hence, 
the single AlAs barrier behaves like a TX filter with a crossover thickness 
Lrx ~ 25A. 

We have also performed analogous calculations for double-barrier struc- 
tures. As mentioned above, a qualitative difference from the single-barrier 
case is that resonant tunneling may happen not only via confined X states 
in the barriers, but also via confined T states in the well. Hence the result- 
ing current-voltage curves display more resonances, both asymmetric Fano 
resonances and "normal" resonances with a Breit-Wigner form. 

Because of the complicated resonant behaviour of Tr and Tx , Rrx is not 
a smoothly increasing function of barrier width. This is already apparent in 
Fig. 6, and in the case of a double barrier the oscillations are even stronger. 
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Even so it is possible to distinguish two regimes, one where mostly Rrx "C 1 
and another where mostly Rrx 3> 1. For the double-barrier structure we 
have evaluated Rrx as a function of AlAs barrier width, for four different 
GaAs well widths: 6, 17, 34 and 56 A. The resulting barrier crossover widths 
are Lrx ~ 14, 22, 25 and 25 A, respectively. This result can be understood as 
follows. When the well is very narrow and consists of only two monolayers of 
GaAs, there is strong coupling between the X states in each barrier. Hence 
the crossover width for one barrier is about half the crossover width in the 
single-barrier case. With increasing well width the coupling between the X 
states in each barrier becomes weaker, and Lrx approaches the value found 
in the single-barrier case. 



VI. Conclusions 

We have presented a framework for the study of tunneling in III-V semi- 
conductor heterostructures. A multiband tight -binding model was used to 
obtain a realistic description of the lowest conduction band of each material. 
We applied the model to GaAs/ AlAs structures where the barrier material 
has an indirect bandgap. Tunneling states ipr and ipx contribute to trans- 
mission. For barriers wider than about ten monolayers of AlAs, tunneling 
via tpx dominates. An applied bias larger than, in the present model, 0.5 V 
enables transfer of T electrons to the X valley in GaAs. The TX transfer 
is stimulated by an indirect-bandgap AlAs barrier. With a conduction- 
band offset of 160 meV between GaAs and AlAs we find a barrier "crossover 
width" Lrx ~ 25 A which separates a regime of negligible TX transfer 
(narrow barriers) from a regime of large TX transfer (wide barriers). 

Finally we would like to comment that the present model is quite gen- 
eral. Given the crystal structure, a set of tight-binding parameters, and 
the conduction-band offsets between the materials involved, one can study 
transport through structures of arbitrary composition. Interesting appli- 
cations, besides the one studied here, could be the polytype type-II het- 



erostructures [37] and the mixing of hole states in a pn junction. One might 
also try to apply the model on the level of the coherent -potential approx- 
imation, first, to enable a study of the effects of alloy scattering in e.g. 
ALrGai-zAs, and second, to see how interface roughness affects the TX 
transfer. 
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APPENDIX A 



From Eq. (2.1) we have 



(nbjk\\\H\n% ! j'k\\) = L~\ j dk x e ik ^ j - j > L/4 (nbk\H\n'b'k), (Al) 



where the matrix elements (nbk\H\n'b'k) are given in table A of Ref. [24]. 
For example, the matrix element between the anion s orbital and the cation 
p z orbital is 

(sak\H\ Pz ck) = V(sa,pc)g 3 {k). (A2) 

Here V(sa,pc) = 4(saR a \H\pcR c ) is one of the, thirteen in all, independent 
TB matrix elements in the basis of symmetrically orthogonalized atomic 
orbitals \nbR), and g% (k) may be regarded as a form factor that is determined 
by the symmetry of the crystal |26[|. One may write 



where d\ = (1,1)ol/4 and di = (1,— 1)ol/4 are vectors in the (y,z) plane. 
Hence ( |A1| ) readily yields 



{sajkv, \H\p z cj'k\\) = V(sa,pc) 



% —* —* % -* -* 

5j+i t j>— sin k^di — Sj^iji— sin k\fa 



(A4) 



for the matrix element in the desired basis. In Eq. ([^) we let U ac represent 
the "hopping" from an anion layer j to a cation layer j — 1, whereas V ac 
represents the hopping in the other direction, from j to j + 1 |]38[| . We then 
have, with s n = ^sinA?||d n and c n = ^ cos k\\d n (n = 1,2), the following 



16 



Compound 


E(s,a) 


E(p, a) 


E(s,c) 


E(p,c) 


E(s*,a) 


E(s*,c) 


V(s,s) 


GaAs 


-8.3431 


1.0414 


-2.6569 


3.6686 


8.5914 


6.7386 


-6.4513 


AlAs 


-7.5273 


0.9833 


-1.1627 


3.5867 


7.4833 


6.7267 


-6.6642 


Compound 


V(x, x) 


V(s,y) 


V(sa,pc) 


V(pa, sc) 


V(s*a,pc) 


V(pa, s*c) 


V(s*,s*) 


GaAs 


1.9546 


5.0779 


4.4800 


5.7839 


4.8422 


4.8077 


0.0000 


AlAs 


1.8780 


4.2919 


5.1106 


5.4965 


4.5216 


4.9950 


0.0000 



Table 1: Empirical matrix elements of the sp 3 s* Hamiltonian in eV for GaAs 
and AlAs. 



elements of our Hamiltonian: 





- E(s,b) 












E(p,b) 










E b = 




E(p,b) 




■ b = c or a , 








E(p, b) 












E(s*,b) \ 








V(s,s)c 2 


-V(sa,pc)c2 


iV(sa,pc)s2 


— iV(sa,pc)s2 







V(pa, sc)c 2 


V(x, X)C2 


-iV(x,y)s 2 


iV(x,y)s 2 


V(pa,s*c)c 2 


u ac = 


—iV(pa, sc)s2 


-iV(s,y)s 2 


V(x,x)c 2 


-V(x,y)c 2 


—iV(pa, s*c)s2 




iV(pa, sc)s2 


iV(x,y)s 2 


-V(x,y)c 2 


V(X,X)C2 


iV(pa, s*c)s2 







-V(s*a,pc)c 2 


iV(s*a,pc)s2 


—iV(s*a,pc)s2 


V(S\S*)C 2 




V(s,s)ct 


V(sa,pc)ci 


iV(sa,pc)si 


iV(sa,pc)si 









—V(pa, sc)c\ 


V(x, x)c\ 


iV(x,y)si 


iV(x,y)si 


-V(pa, s*c)ci 




Vac = 


—iV(pa, sc)si 


iV(x,y)si 


V(x, x)c\ 


V(x,y)ci 


iV(pa, s*c)si 






—iV(pa, sc)si 


iV(x,y)si 


V(x,y)ci 


V(x,x)ci 


iV(pa, s*c)si 









V(s*a,pc)ci 


iV(s*a,pc)si 


iV(s*a,pc)s\ 


V(s*,s*) Cl 





In addition U r . 



U\ c and V a 



, (A5) 

Vj„. Numerical values for GaAs and AlAs, 



taken from Ref. [24], are collected in Table 1. 



APPENDIX B 

In this Appendix we will follow a general approach and obtain the 
electric-current density in terms of nonequilibrium Green functions. Using 
socalled "surface Green functions" we will show how to make contact with 
the transmission-probability approach presented in Sec. IV. The present 
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treatment is closely related to the one introduced by Caroli et al. 3£]. For 



a system of noninteracting electrons, the total current density J is given 
in terms of the single-particle current-density operator J as a many-body 
quantum average over occupied single-particle states, and a statistical aver- 
age over all possible many-body configurations. Thus, 

J = E^HEft»W(^A), (si) 

{m} s 

where P(m) is the probability of having the many-body configuration {m}, 
and p m (s) is the probability that the single-particle state \4> s ) is occupied in 
the given many-body configuration \p m (s) equals or 1]. Since the reduced 
density matrix p is defined as 

p=£ P MEft»Wl60fa'l> (B2) 

{m} s 

we may also write 

J = tr(pJ), (B3) 

where tr denotes the trace operation. The connection to Green functions is 
now transparent since 

P = f^ G< ( E )> (B4) 



where we follow the notation in e.g. Ref. ||40(| . The correlation function 
G < and its hole "counterpart" G > are related to the retarded and advanced 
Green functions G R and G A via 4C] 



G < + G > = i(G — G ). (B5) 

In general only two of the four functions in ( |B5| ) are independent since one 
also has the connecting identity 



G 



R 



* . (B6) 



Eqs. ( |B3| ) - (B6) constitute the proper starting point for evaluation of the 
electric-current density in a system of noninteracting electrons. 

An operator J for the electric-current density is found by comparing the 
continuity equation and the Heisenberg equation of motion. We also use 



(B4) and obtain the result 



Jm,m+1 — fi j 2lT l (27r) 2 ^ r ^"* +1 ' mltm ' m+1 G m,m+l u m+l,m) (B7) 
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for the net current density flowing between layers m and m + 1 (in the 
positive x direction). In ( |B7| ) the matrix elements of G < depend, of course, 
on E and A;m, and the TB "hopping" elements u mjm ±i depend on ku in a way 
which is determined by the crystal structure (see Appendix A). Eq. ( |B7| ) is 
valid for systems that have translational invariance in the parallel plane and 
are characterized by a nearest-neighbour TB Hamiltonian. 

In order to connect the Green-function approach and the transmission- 
probability approach, we must show that 



ztr(G, 



< 



m+l,m u m,m+l 



T aP ifFD(E - Hl) ~ fFD{E ~ HL + < 

a/3 

(B8) 

when we have a structure connected to semiinfinite perfect leads in equilib- 
rium at chemical potentials hl and hr = //£, — eV. The following ingredients 
are required to rewrite the left-hand side of this equation: 

• The equilibrium expressions for G < and G > in a homogeneous system 
at chemical potential fi plj]: 



9 < (E) 
9 > (E) 



g R (E)-g A (E) 
g R (E)-g A (E) 



Jfd(E - fx), 
[l-f FD (E-ri] 



(B9) 



(We use lowercase letters for Green functions that describe a system 
in equilibrium.) 



The rules, derived by Langreth and Wilkins [42], for handling products 
of two or more Green functions: 



(AB) R 
(AB) A 
(AB)< 
(AB)> 



A R B R , 
A A B A , 

A R B<+A<B A , 
A R B> +A>B A . 



(BIO) 



• The Dyson equation and the recursive Green-function technique, de- 
scribed in detail in the Appendix of Ref. [3C|. 

• The unit matrices of Eq. ( |3.11 ). 

Straightforward algebra then yields 

(a c |C/ ca 2Im7+C/ ac | / a c ){/9|G^ +li0 |//}(M o |C/ ac 2Im 7 -C/ co |^) x 

a(3fiv 



{v\G A N+1 \a}[f FD (E - n L ) - f FD (E - m + eV)] 



(Bll) 
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for the left-hand side of flB8|), where we have chosen to let the cation and 
anion layers of "site" N + 1 represent m and m + 1, respectively. Here 
subscripts c or a on the states a, f3, /i and v denote the cation or anion 
parts, respectively. Further, 2Im7 = i{^ A — j R ) |f43| , and 7^(7^) is the 
surface Green function, i.e., the local element at the surface layer of the 
Green function for a semiinfinite crystal extending from layer N + 1, a (0, c) 
to +00 (—00). Remember, sites N + 1 and are defined as the first sites 
of the semiinfinite crystals at constant potential energy on each side of the 
heterostructure. From the recurrent relations for the surface Green functions 
[p0[ it is clear that 7+ is the same on site N + 1 as in +00, and 7" the same 



on site as in —00. In other words, 7+ and 7 C in ( Bll ) reflect asymptotic 
properties of the left and right lead, respectively, which allows for the use 
of © in deriving ( [Bill ). 



In ( |B11 ) only a few terms contribute to the sum. First, there can only 



be a contribution from extended states with real values of k. Clearly, the 
result of ( Bll| ) must be unchanged if we choose sites, other than N + 1 



and 0, further into the asymptotic regions. However, evanescent states (5 
or [i would give a modulus that decreases exponentially with 

increasing i or decreasing j (i > N + 1, j < 0). Second, only states with 
positive velocity contribute to the sum since the retarded Green function 
propagates a scattering state forward in time. Finally, since G A = [G R y , 
the states v and a must also be extended states with positive velocities. 

What remains is to evaluate the matrix elements (a c \U ca 2Im'y^U ac \(5 c ) 
and 

(^ a |[/ ac 2Im7~[/ ca |z/ a ). This task is accomplished by rewriting matrix el- 
ements of the velocity operator v [see Eq. (|3.5|) 1 with help of the Dyson 
equation and the recursive Green-function technique. The result is 



(a c |C/ ca 2Im 7 +C/ ac |/? c ) = -%Ra\v\(3) = -^v a 5 af} (v a > 0), 
(fi a \U ac 2lm-f-U ca \u a ) = ^{^\v\v) = HT^V < 0). 



(B12) 

Here the relation (a|-0|/?) = v a 5 a p follows from current conservation. Al- 
though states with negative velocity do not contribute to the sum in ( |B11 ) , 



there corresponds to each state \x (with < 0) a "time-reversed" state \j! 
(with = — Vfj, > 0) for which one has 

Ah 

(^|C/ ac 2Im 7 -C/ ca |^ a ) = V . (B13) 
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Finally, collecting all our knowledge, the sum in ( |B11| ) may be written as 



V —v a S a /3{P\G§ +10 \a}—v fM 5, xu {i'\GQ N+1 \a}[f FD (E - fi L ) - Jfd(E - hl + eV)] 
aPw aL aL 

E f — ) v«vp\{(3\G% +l . \a}\ 2 [f FD (E - n L ) - f FD {E - y, L + eV)] 



= E T «p[fFD(E - n L ) - f FD (E - ii L + eV)] , 

a/3 



(B14) 



which demonstrates the validity of Eq. (B8), and hence the equivalence of 
the transmission-probability approach in Sec. |y| and the Green-function 
approach. 



* Present address: Institutt for fysikk, Norges Tekniske H0gskole, Univer- 
sitetet i Trondheim, N-7034 TVondheim-NTH, Norway. 
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Figure Captions 



FIG.l. Projection of the zincblende lattice on the (100) direction. Solid and 
open circles represent cation and anion layers, respectively. The interlayer 
distance is ol/4, i.e., one forth the lattice constant. Indicated are also 
the "on-layer" energies E c and E a , and the interlayer hopping matrices 
U ca , U ac , V ca and V ac . 

FIG. 2. Parts of the complex band structure for GaAs and AlAs. The solid 
lines are the lowest conduction band of the two materials. For AlAs we have 
also included the imaginary parts, Kr and kx, of the analytical continuation 
below the minima at the T and X points, respectively. [To be precise, the 
conduction-band minimum of AlAs is not exactly at the X point in the 
present model, but slightly below, at k ~ 0.85 x (27r/aL).] Zero energy is 
taken to be at the top of the valence band in GaAs, and with the present 
model and parameters GaAs has an energy gap of 1.55 eV. 

FIG. 3. Spatial decay into AlAs of the wave functions ipr (thick lines) and 
ipx (thin lines) at energies 10 meV (solid lines) and 100 meV (dashed lines) 
above the GaAs conduction-band edge £j(GaAs). The distance from the 
GaAs/ AlAs interface is given in terms of the number of AlAs monolayers 
N^. The decay lengths corresponding to these curves are Kp 1 = 0.95ol 
and 0.98<2l ; k^ 1 = 0.92ol and 1.78cil- Here the first (second) number 
refers to the solid (dashed) lines, and = 5.66 A is the lattice constant in 
AlAs. The inset illustrates the V point conduction-band minimum which 
represents a potential step between iVj, = and iV& = 1. The wave functions 
are normalized to an incoming plane wave of unit amplitude (see Sec. III). 

FIG. 4. Tunneling of an electron through a single AlAs barrier with applied 
bias V. The profiles of the T and X point conduction-band minima are 
represented by the solid and dashed lines, respectively. The lowest conduc- 
tion band in GaAs [in the (100) direction] is sketched on both sides of the 
barrier, and the dotted portions indicate initially occupied states on the left 
side and available transmission states on the right side. Shaded areas denote 
the filled equilibrium Fermi sea in the emitter and collector contact. The 
transmitted electron may contribute to the "direct" current density Jr or 
the "transfered" one Jx ■ The tunneling process is described in detail in the 
text. Relaxation processes that bring the transmitted electrons to thermal 
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equilibrium are indicated with scattering rates r rr and t xt . 

FIG. 5. Direct and transfered current-density contributions (thick lines), 
Jr and Jx, for a single AlAs barrier of width 40 A. The corresponding 
transmission coefficients (thin lines), Tr and Tx, at E = Ep = 10 meV 
and fey = reflect to a good approximation both the qualitative behaviour 
of Jr and Jx, and also their relative contributions to the total current 
density. The left vertical axis represents Jr and Jx in A/cm 2 ; the right axis 
represents Tr and Tx which are of the order of 10 -4 in this case. 

FIG. 6. Rrx, as defined in Eq. (5.2), for a single AlAs barrier as function of 
barrier width, the latter given in terms of Nf,, the number of monolayers of 
AlAs. The solid line connecting the data points is nothing more than a guide 
to the eye. Rrx gives a quantitative measure of the relative contribution 
of Jx and Jr to the total current density: Rrx *C 1 implies negligible TX 
transfer; Rrx 3> 1 implies large TX transfer. 
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